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Abstract 

We develop a technique to study relativistic perturbations in the generahsed brane cosmological 
scenario, which is a generahsation of the multi-fluid cosmological perturbations to brane cosmology. 
The novelty of the technique lies in the inclusion of a radiative bulk which is responsible for bulk- 
brane energy exchange, and in turn, modifies the standard perturbative analysis to a great extent. 
The analysis involves a geometric fluid - called the Weyl fluid ~ whose nature and role have been 
studied extensively both for the empty bulk and the radiative bulk scenario. Subsequently, we find 
that this Weyl fluid can be a possible geometric candidate for dark matter in this generalised brane 
cosmological framework. 
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I. INTRODUCTION 



During the last few years braneworld gravity has emerged as a more general theory of 
gravity, mainly due to the possibility of explaining the gravitational phenomena observed in 
the four dimensional universe from a broader perspective [l, 2]. Subsequent developments 
of the theory in the cosmological sector came as an inevitable outcome since the chal- 
lenges any theory of cosmology, be it a theory based on General Relativity or any other 
phenomenologically motivated theory, faces in explaining predictions from the highly accu- 
rate observational data 0, O, gI- In spite of great complications involved, the cosmological 
aspects of this scenario did show some promising features. To mention a few, brane cosmol- 

n 

ogy naturally gives rise to singularity-free bouncing and cyclic universes [7|. Also, in this 
theory, the universe does not need any special initial condition for the inflation to start so 
that the isotropy is built in the theory ^]. Even the possibility of inflation without any 4D 
inflaton field is in vogue [oj] . Brane cosmology thus results in interesting physics which needs 
further investigations. 

In this scenario, the bulk spacetime is either AdSs flO^ or a generalised version of it. The 
generalised global structure depends upon whether the bulk has only a cosmological constant 
or there is any non-standard model fields minimally or non-minimally coupled to gravity or 
to brane matter. When the bulk is empty consisting only of a cosmological constant, the bulk 
metric in which an FRW brane can be consistently embedded, is given by a 5-dimensional 
Schwarzschild-Anti de Sitter (Sch-AdSs) or a Reissner-Nordstrom Anti de Sitter (RNAdSs) 



black hole 



12] . A subsequent generalisation of this scenario can be obtained when 



d, resulting in a Vaidya- 
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the bulk is not necessarily empty but it consists of a radiative fie 
Anti de Sitter (VAdSs) black hole for the bulk metric lal,!!^, 14, 3, 
A 'black hole in the bulk' scenario provides us with a novel way of visualising cosmological 
phenomena on the 4D universe. In this scenario, the brane is moving in the bulk, with its 
radial trajectory being identified with the scale factor of the 4D world, so that the expansion 
of the universe is a realisation of the radial trajectory of the brane in the bulk. 

The most notable contribution from bulk geometry on the brane is, perhaps, an additional 
term in the Friedmann equations, which arise from the projection of the bulk Weyl tensor 
onto the brane. The precise role of this term, compatible to FRW background on the brane, 
is to supply a geometric perfect fluid whose nature is governed by the contents of the bulk (in 
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turn, bulk geometry) we choose. For an empty bulk, it is radiation-like and is called the dark 
radiation. There is extensive study in the literature either by setting it to zero for practical 
purpose or by attributing a very small value to it, constrained by Nucleosynthesis data 



[< 3% of total radiation energy density of the universe) 0, |23|]. Examples include metric 



)ased perturbations [2J], density perturbations on large sea 



25 1, or including its effect s |26[| . curvature perturbations |27| and the Sachs- Wolfe effect 



es neglecting dark radiation 



281], vector perturbations [23], tensor perturbations 30| and CMB anisotropics [31]. In 



all the cases, the effect has been found to be modify the standard analysis very little, as 
expected from its radiation-like behaviour. 

On contrary, when the bulk is not necessarily empty, the nature of this entity is no longer 



radiation-like, rather it depends 



VAdSs bulk geometry 
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le contents in the bu^ 
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k, which is reflected by the 



33|, 



341 ]. It is thus important 



to determine the nature as well as the role of this entity, called in general the Weyl fluid, in 
the cosmological dynamics and perturbations, and find if this scenario has some advantages 
over others. A recent work 2l| has shown its significance as a possible dark matter candidate 
)y Newtonian analysis of perturbations, followed by some confrontation with observations 



351 ] . However, as in GR, the Newtonian analysis of gravitational instability is limited in the 



sense that it fails to account for the perturbations on scales larger than the Hubble radius. 
One needs relativistic analysis valid at super-Hubble scales as well. Further, in order to 
test braneworld scenario observationally, we need a complete description of the evolution 
of density perturbations in the most general brane cosmological scenario provided by this 
VAdSs bulk. With these motivations, we develop here a technique for relativistic density 
perturbations valid for this generalised brane cosmology, which will act as a natural extension 
of the covariant perturbations of General Relativistic framework [sgI to braneworld scenario. 
We further show in the subsequent discussions that the Weyl fluid can play a crucial role 
in late time cosmologies as a geometric candidate for dark matter albeit its actual material 
existence. 
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II. BRANE DYNAMICS WITH WEYL FLUID 



As mentioned, we shall concentrate on the most general bulk scenario, for which the bulk 
geometry is given by a Vaidya-anti de Sitter metric 



-f{r, v) dv^ + 2dr dv + r^dEj 



(2.1) 



where S3 is the 3-space having flat, spherical or hyperboloidal symmetry, /(r, v) 



k - 



A5^2 
6 



m,(v) 
r2 1 



^r" — and m{v) is the resultant of the variable mass of the Vaidya black hole and 



radiation field. This t ype of bulk can exchange energy with the brane as a null flow along 



the radial direction 
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3, 3: 
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22l |. Consequently, the brane matter 



conservation equation is modified to 



p+ 3-(p + p) = -2^ 

CI 



(2.2) 



where ip is the null flow characterising the VAdSs bulk by the radiation field of a null dust 
'^AB^ — i^lAQs, which leads to the above equation by using 



.0 rbulk „A B 



(2.3) 



(where n-^ are the normals to the surface), which gives 



-2V'M, 



(2.4) 



(m^ are the unit velocity vectors), and readily leads to Eq (12.21) . This modified conservation 



equation, with the help of the Bianchi identity on the brane V^G^^, 
constraint equation 



0, leads to another 



A ^ '''^'^ + 3 



a 



u, 



(2.5) 



where A is the brane-tension and S^u and S^i, are, respectively, the projected bulk Weyl tensor 
and the quadratic contribution from the brane energy-momentum tensor to the Einstein 
equation on the brane 2l|. The above equation governs the evolution of the Weyl fluid p* 
(so named since it is a fluid-like contribution from the bulk Weyl tensor to the brane). For 
FRW geometry on the brane, this is given by 



p* + 4-p*=2^- 
a 6 



2 /K5 



a 



(2.6) 



so that this quantity evolves as ly, |21 1 



P 



CO 



1 



oc 



(2.7) 



which gives a general, physically relevant behaviour for the Weyl fluid. Here, r is the proper 
time on the brane. Obviously, contrary to the Sch-AdSs bulk, here C(r) is evolving, and 
consequently, the Weyl fluid no longer behaves like radiation. To a brane-based observer, 
the cosmological dynamics is now governed by an effective perfect fluid, the components of 



which are given by 



15|, 



2l| 



eff C{t) 



(2.8) 
(2.9) 



The anisotropic components of the Weyl fluid, viz., g* and vr*^^ vanish, in order that the 
VAd^s bulk be compatible to FRW geometry on the brane. The Friedmann equation and 
the covariant _^aychaudhuri equation, expressed in terms of these effective quantities, are 
respectively 



22 



^2 ^ ^^cff + ^ _ 



(2.10) 
(2.11) 
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In the brane-based Newtonian analysis of perturbations by considering small fluctuations 
of the effective density p'^^{lt, r) = p'^^{t){1 + 6'^^{lt, r)) and the effective gravitational 



potential $ 



eff/ 



X , r 



$Q + 4> on the hydrodynamic equations for this effective perfect 



fluid, one obtains for a barotropic fluid a single second order equation in terms of Fourier 
mode 



+ 2 



a dS 



eff 



a dr 



sreff 







(2.12) 



where cf^^ is the square of the effective sound speed [2], [21I . The above perturbation equation 
of the effective fluid can account for the required amount of gravitational instability if the 
Weyl density redshifts more slowly than baryonic matter density, so that it can dominate 
over matter at late times, which is realised when 1 < a < 4 in Eq (12. 7p . Now, for late time 
behaviour, we can drop the quadratic terms in equations (12. 8p and (12. 9p so that the effective 
density is given by p'^^ = p*-^-* + p* which is now constituted of the usual matter (baryonic) 
density p^''^ and an additional density contribution from the Weyl fluid. This Weyl density. 



being geometric, is essentially non-baryonic. Consequently, we can decompose Eq (12.120 to 
get the individual evolution equations for the perturbation for each of the fluids 

^ + = 47rGpW5(^) + A^Gp*5* (2.13) 

d^ ^ ^ad^ ^ ^ 47rGp(^)5W (2.14) 

dr^ a dr 

where d'^^^ and 5* are the fluctuations of baryonic matter and Weyl fluid respectively. With 
-C the relevant growing mode solutions are given by 21 1 

5*{z) = S*{0){l + z)-' (2.15) 
S('\z) = 6*{z) (l - ^i^) (2.16) 



1 + Zn. 

with the input that the late time behaviour of the expansion of the universe in RS II is the 



same as the standard cosmological solution for the scale factor 1^, 32| where the scale factor 
is related to the redshift function by a oc (1 + 2)^^. 

The solutions reveal that at a redshift close to zn, the baryonic fluctuation 5^^^ almost 
vanishes but the Weyl fluctuation S* still remains finite. So, even if the baryonic fluctuation 
is very small at a redshift of zn ~ 1000, as confirmed by CMB data [4], the fluctuations 
of the Weyl fluid still had a finite amplitude during that time, whereas at a redshift much 
less than z^ the baryonic matter fluctuations are of equal amplitude as the Weyl fluid 
fluctuations. This is precisely what is required to explain the formation of structures we see 
today. Thus, the Newtonian analysis of perturbations on the brane is capable of explaining 
structure formation (within its limit) by Weyl fluid, devoid of any material existence of dark 
matter. Hence, the Weyl fluid acts as a possible geometric candidate for dark matter. 



III. RELATIVISTIC PERTURBATIONS WITH WEYL FLUID 

The Newtonian analysis depicted so far turn out to be an useful tool to study perturba- 
tions on the brane after the universe enters the Hubble length. A more complete picture 
can be obtained if one studies relativistic analysis of perturbations, which include the evo- 
lution of the universe at the super-Hubble scale as well. In this section we shall develop a 
multi-fluid perturbative technique in order to discuss relativistic perturbation relevant for 
brane cosmology. This will be carried through in the subsequent sections for the purpose of 
analysis for different braneworld scenarios. Our basic motivation in the attempt to develop 
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a multi-fluid perturbative technique is governed by the reahsation obtained from Newtonian 
analysis that, contrary to the Sch-AdSs bulk scenario, the Weyl fluid may not be that much 
insigniflcant so as to neglect its effects at late time, if we have a general bulk geometry. 
Consequently, in a general brane cosmological scenario, along with baryonic matter, the 
universe consists of a considerable amount of Weyl fluid as well. 

Before going into the details, let us jot down here the major points in addressing rela- 
tivistic perturbations on the brane. 

• Here the cosmological dynamics is governed by a two-fluid system. One of the compo- 
nents of the system is a material fluid p*^''^ - the baryonic matter content on the brane. 
The second component is a geometric fluid p* - the Weyl fluid. The total (effective) 
density for the system on the brane is given by p^ff = pW +p*. 

• Though there are two components of the effective fluid, the Weyl fluid being a geo- 
metric entity, there is a single material fluid in the analysis. As a results, there will 
be no entropy perturbation as such. 

• For the same reason, there is no peculiar velocity for the Weyl component, leading to 



• The anisotropic components of the Weyl fluid being absent so as to flt it into an FRW 
background, we will set = = vr*^ right from the beginning. As a result, each 
component of the two-fluid system behaves individually like a perfect fluid, resulting 
in a perfect fluid behaviour for the effective fluid as a whole. 

• These two fluids interact and exchange energy between them, which is governed via 
the bulk-brane energy exchange and the backreaction of the system on the brane. 

• Since there is energy exchange between these two fluids, the conservation equation for 
each individual is now modifled. These modifled forms of the conservation equations 
have been explained in equations (12.21) and (12.61) . 

Because of the interaction between the two fluid components, each of the two modifled 
conservation equations can be written in terms of the contribution from the interaction as 



V* = 0. 




(3.1) 
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where 9 = 3^ is the volume expansion rate, a superscript (?) denotes the quantities for the 
2-th fluid and /'■*•' is the corresponding interaction term. It readily follows from equations 
fl2.2p and fl2.6p that the interaction terms, when written explicitly, are given by 

= (3.2) 

r = 2^l:-^(^y (i^ + 3^ii] (3.3) 



3 \ K J \ a 

For relativistic perturbations, we express the densities of each of the contributing fluid 
components in terms of dimensionless parameters as 

with the first one for baryonic matter and the second one for Weyl fluid. Considering 
the nature of the Weyl fluid as discussed in the previous section, we find that the density 
parameter for the Weyl fluid is given by 

^P* = (3.5) 

where Cq is the onset value for the Weyl parameter C(r). 

For completion, we mention here that there can, in principle, appear two more dimen- 
sionless parameters, one each for the cosmological constant and the brane tension arising in 
the brane cosmological context. They are 



with the total density satisfying the critical value 

^ot = J2^i = ^pC) +np, + nA + nx = i (3.7) 

i 

Here, and throughout the rest of this article, we have considered a spatially flat universe 
with k = 0. In Eq (13. 6p . the first one is relevant if one considers cosmological constant in 
this brane universe while studying the expansion history of the universe whereas the one 
due to the brane tension is relevant in the high energy early universe (inflationary) phase 
but is negligible for low energy late time phenomena such as structure formation. Thus 
the baryonic density and the Weyl density are the only two relevant contributions in the 
scenario being discussed here. In what follows we shall restrict ourselves to the discussion 
of the Einstein-de Sitter brane universe for which = leading to f2tot ~ + ^p* = 1. 
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We now express the comoving fractional gradients of the effective density and expansion 
relevant in the brane cosmology as 

^ = ^D,p^' (3-8) 
Zf = aD,Q (3.9) 

< = ^D,p^' (3.10) 

As already discussed, both baryonic matter and Weyl fluid behave individually as perfect 
fluid components, which means the effective flux arising from the peculiar velocities of each 
component vanish to zero order, confirming that the perturbations considered here are gauge- 
invariant at the first order. 

With the above notations, the linearised evolution equation for the density perturba- 
tions in the braneworld is obtained by taking spatial gradient of the modified conservation 
equations. After linearisation, it turns out to be 

y pM J I" ^ ' pW(i + iyeff) pW p(i) ^ ' 

where ti'^*-' = p*^*^/p*^*^ is the equation of state for i-th fluid and c^*-*-* = p'-^-'/p''*'' is the sound 
speed squared for that species, with the corresponding quantities for the effective total fluid 
are, respectively, 

1 



w 



P 



c 



Y,cT^P^'\^ + w^'n (3.13) 



In the relativistic perturbations, contrary to the Newtonian analysis, we further have an 
evolution equation for the effective expansion gradient, which depends on the effective fluid. 
This is obtained by taking spatial gradient of the modified Raychaudhuri equation ( 12. lip 
and is given in the braneworld scenario by 

zf + 2HZf = -f /"A- - ^D.D^Af + j^c^Af - a.lD,,p (3.14) 

It should be mentioned here that since the evolution of the expansion gradient is depen- 
dent on the curvature perturbations, the later should not remain strictly constant in this 
multi-fluid perturbation scenario. However, though there is a significant energy exchange 
between brane matter and Weyl fluid at early times, we shall see from the next section that 



the energy exchange between the two fluids are almost in equilibrium at late times, so that 
the local curvature perturbations can safely be considered to be constant for all practical 
purpose. One should, however, consider the variation of this term while analysing infla- 
tionary phase for an instance. We follow this argument right from here in order to avoid 
mathematical complexity. 

As in GR, we find that while discussing perturbations in brane cosmology, it is advan- 
tageous to express the above equations in terms of covariant quantities. These density 
perturbations are governed by the fluctuation of the following covariant projections 

= a D^'Af (3.15) 
A"^ = a D^Af (3.16) 
Z^^ = a VZf (3.17) 

Consequently, the covariant density perturbation equation and expansion gradient on the 
brane, when expressed in terms of the above covariant quantities, are obtained straightaway 
from equations (13.111) and (13.141) . They are given by 



A« = 3W) - — A« - (1 + u;«)Z^*^ - 



.eff 



_^^effy 



+ -TsD^I^^ 3.18 

^eff + 2i7Z^ff = --p^'^A'^f^ - ^^D^A^'^ + ^^^cf ^A^f^ - aVD^^ (3.19) 

In deriving the above covariant perturbation equations, we have considered those kind of 
perturbations for which, like the unperturbed Weyl fluid, the anisotropic stresses and fluxes 
for the perturbed Weyl fluid are also vanishing. 

These set of equations provide the key information about the perturbation in brane 
cosmology. In the subsequent section, we shall try to analyse these relativistic perturbation 
equations and obtain possible consequences. 



IV. SOLUTIONS AND ANALYSIS 

A. Empty bulk : Non-interacting fluids 



Let us now discuss the special scenario when the bulk is empty for which the VAdSs 
bulk reduces to Sch-AdSs. In this case, there is no question of energy exchange between the 
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brane and the bulk. Consequently, there is no interaction between brane matter and Weyl 
fluid as such, which reveals from Eq (13. ip the fact that the individual conservation equation 
for each of the components are preserved. Thus, the Weyl fluid evolves in this case as 

p* oc a"^ (4.1) 

with the Weyl parameter a now being zero, so that for empty bulk, the Weyl fluid behaves 
like radiation, for which this is called dark radiation. 

Since in this case, there is no interaction between the two components of the effective 
fluid and also, there is no null flow from the bulk to the brane (or vice versa), we can drop 
the interaction terms and the terms involving ip in the analysis. As a result, the covariant 
perturbation equations (13.181) and (I3.19P are vastly simplified. They are now given by 

A» = 3i7w(^)A« - (1 + ti;W)Z^^ (4.2) 
Z-ff + = --p^f^A<=fif - ^^D^A-ff (4.3) 

Taking the time derivative of the above two equations and combining them, we obtain 
the evolution equations for density perturbations of the two fluids 

AW + 2i7A(^) = ^p'^A""^ (4.4) 

4k2 / k'^\ 
A* + 2HA* = p'^A"^ + A* 21/2 + hA* (4.5) 

3 2 y ^ / 

where the first equation is for baryonic matter while the second one for dark radiation. 

Recall that the amount of dark radiation is constrained by the Nucleosynthesis data to 
be at most 3% of the total radiation density of the universe. So, it redshifts at a faster rate 
than ordinary matter on the brane so that the matter on the brane becomes dominant on 
the Weyl fluid at late time. Hence, it is expected that the dark radiation does not play 
any significant role in late time cosmologies. It is obvious from the fact that in this case, 
^ ^l*, which when put back into the above equations, leads to A^^^^ ^ A*, so that the 
dark radiation fluctuation does not contribute substantially at late times. The Sch-AdSs 
bulk scenario thus fails to explain structure formation with only baryonic matter and dark 
radiation. One needs cold dark matter in the theory and the dark radiation can, at best, 
slightly modify the standard perturbative analysis. 
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B. General bulk : Interacting Weyl fluid 



The general scenario, however, is different from the empty bulk case since now the Weyl 
fluid exchanges energy with brane matter through interactions and is the dominant contribu- 
tion of the effective fluid in the perturbation equations. Here, the evolution of perturbations 
for the individual fluids are governed by equations fl3.18p and (13.191) . which now include the 
effects of the interaction terms as well as of the effect of null radiation through the term 
involving ip. With these inclusions, the equations become a bit too complicated and it is al- 
most impossible to have an analytical solution from these complicated equations. However, 
the equations turn out to be tractable if we incorporate certain simplifications following 
physical arguments, without losing any essential information as such. The simplifications 
we incorporate are as follows: 

• The null flow from the brane to the bulk ip is a, function of time only. This means 
that we are considering only the time-evolution for the null radiation, at least on the 
brane, which is relevant for its late time behaviour in perturbation analysis. 

• The energy exchange between the two fluids is in equilibrium, i.e., the energy received 
by the Weyl fluid is the same as the energy released by brane matter, so that = 
j{b) _|_ /* = Q Hence, no extra energy is leaked to the bulk from the brane at late 
time (though at early time there may be some leakage of energy from the brane to 
the bulk). This basically describes the late time behaviour, consistent with the fact 
that the standard evolution history (scale factor) are regained in this scenario at the 



"matter-dominated" era 
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32|. 



Now, we have shown that in this generalised braneworld scenario, the Weyl fluid, in 
general, evolves as (Ref Eq (12.71) ) 

p* = Coa-('-") (4.6) 

with the parameter a in the range 1 < a < 4 so that it is the dominant contribution in 
the two-fluid system. The energy exchange between the components of the system being in 
equilibrium, we find from Eq (12. 7p that the Weyl fluid now behaves as 

p* oc (4.7) 

with the parameter a = |. This readily suggests that the Weyl fluid actually redshifts more 
slowly than ordinary matter and hence, can dominate over matter at late times, reflecting one 
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of the fundamental properties of dark matter. This also provides a more stringent bound for 
the value of a from theoretical ground alone (which was predicted from Newtonian analysis 
to fall within 1 to 4). 

We now take the time derivative of the covariant perturbation equations fl3.18p and (13.191) . 
and rearrange terms so as to obtain a single second order differential equation for each of 
the fluids. Thus, the equation describing evolution of scalar perturbations of matter on the 
brane turn out to be 



(4. 



whereas the scalar perturbation equation for the Weyl fluid on the brane is given by 

3 2 1 + w^f^ V P* 3 p* } 

- + A- [lE^ -±JJ1±-^±\^K'{h-^±-'^^\ (4.9) 

l^w"^ p* \ 2 p* p* j \ p* 3 J 

Recall from the discussions following Eq ( 14. 7p that in this scenario the Weyl fluid is the 
dominant component of the effective fluid. Consequently, the evolution equation for the 
Weyl fluid at late times is radically simplified by using << A* since the Weyl fluid is 
now the dominant contribution. With the energy exchange between the two fluids being in 
equilibrium, the expression for the null flow further simplifies the above equation so that it 
can now be recast in the following form 

A* + -A* - f- + ^) A* = (4.10) 



t \t t\ 

where the constants A, B, C are readily determined from the constraint equations. These 
constants are given by 

A-um\a^T (4.11) 



3 2\p*J\3Ho 



The above equation (14.101) for A* turns out to be somewhat tractable. One of its solutions 
is given by 



A* ~ t2" -Bessel/ 



Vl-2A + A'^ + 4C, 2VBVt\ (4.14) 
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The above solution, consisting of a Bessel function, is found to be a growing function. 
Therefore, the evolution equation for the Weyl fluid, indeed, shows a growing mode solution, 
which is required to explain the growth of perturbations at late times. Thus, the relativistic 
perturbation theory relevant in brane cosmology gives rise to a fluid which is very different 
from dark matter in origin and nature but has the potentiality to play the role of dark 
matter in cosmological context. It is worthwhile to note that, to a brane-based observer, the 
nature of the Weyl fluid is determined from bulk geometry arising from the radiation flow 
in the bulk. That is why the Weyl fluid can be treated as a geometric candidate for dark 
matter. 

The following figure depicts a qualitative behaviour of the growth of Weyl fluid pertur- 
bations with time. The figure once again shows that the evolution of perturbations of Weyl 
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FIG. 1: Growth of Weyl fluid perturbations with time 



fluid is very different from cold dark matter (CDM), which makes the theory distinct from 
standard analysis with CDM. We, however, note that since the dynamics here is completely 
different from the standard one involving CDM, one cannot comment conclusively on the 
merits/demerits of Weyl fluid over CDM right from here. One has to reformulate and esti- 
mate different cosmological parameters in this context and confront them with observations 
for a more conclusive remark. For example, the relation of the transfer function with the 
potential will now be replaced by a novel relation with the effective potential discussed in the 
brane cosmological context [21,]. As a result, the variation of the growth function with the 
scale factor may not be the same as usually needed in the standard cosmological paradigm. 
It is to be seen if this analysis of perturbation with the Weyl fluid fits in this new, brane 
cosmological framework, which is not a trivial exercise, we suppose. The interested reader 
may further refer to 37| for an overall view on how different cosmological parameters are 
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developed in a specific theoretical framework. 

However, even at this stage, our model does show some agreement with observational 



results. From the recent studies on confronting braneworld models with observations [35 1 
by obtaining the luminosity distance for FRW branes with the Weyl fluid, it is found that 
a certain amount of Weyl fluid with 2 < a < 3 is in nice agreement with Supernovae data. 
From the relativistic perturbations discussed in this article we have found a specific value for 
a, namely a = |, which falls within this region. Thus the braneworld model of perturbations 
fits well in this observational scenario. We hope an extensive study in this direction will 
lead to more interesting results to make a more conclusive remark. 



V. SUMMARY AND OPEN ISSUES 

In this article, we have developed a technique for relativistic perturbations valid for a gen- 
eral brane cosmological scenario. The essential distinction of our analysis from the studies 
on brane cosmological perturbations available in the literature is that, here the geometri- 
cal effect of the bulk on the brane - the so-called Weyl fluid - plays a very crucial role in 
determining the nature of the evolution of density perturbations. This is materialised from 
the realisation that in the general brane cosmological scenario obtained from Vaidya-anti 
de Sitter bulk, the Weyl fluid plays a significant role in controlling the dynamics on the 
brane, contrary to the earlier results based on dark radiation. Our results are, in a sense, 
a generalisation of the multi-fluid covariant perturbation formalism in brane cosmological 
framework. Further, we have solved the perturbation equations and found that the pertur- 
bation of the Weyl fluid grows at late time, and thus, this component of braneworld gravity 
plays a significant role in late time cosmology to act as a possible geometric candidate for 
dark matter. We have discussed some of the implications of fluctuations involving it and 
have mentioned some observable sides of this model as well. 

An important issue is to fit this theoretical model with current observational data. Re- 



cently there has been some progress in this direction [35|] . An extensive study on confronting 
this braneworld model with observations in a more rigorous method can provide us with nec- 
essary information on the merits and demerits of the formalism. To this end, a thorough 
study of different parameters related to cosmological perturbation is to be performed. As 
mentioned in the previous section, the different cosmological parameters need to be refor- 
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mulated in this framework. The next step is to estimate them and confront them with 
observations. For example, it is to be seen if the power spectrum, redefined in this paradigm 
with the Weyl fluid acting as a dark matter candidate, fits with the highly accurate obser- 
vational data. 

Further, analysis of different types of metric-based perturbations, namely, scalar, vector 
and tensor as well as related issues like CMB anisotropy, Sachs- Wolfe effect etc has to be 
be studied in details in this brane cosmological framework with a significant Weyl fluid. An 
extensive study in this direction is essential, which we hope to address in near future. Also, 



to apply this formalism in the framework of branworld models of dark energy 38|J remains 
as another interesting issue. 
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